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EXTENDED ACTIVE DISTURBANCE 
REJECTION CONTROLLER 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

This application claims the benefit of U.S. Application No. 
60/965,359 filed 20 Aug. 2007. This application is a continu- 
ation in part and claims the benefit of the following Applica- 
tions: U.S. application Ser. No. 12/067,141 filed 17 Mar. 
2008, and issued as U.S. Pat. No. 8,060,340 on 10 Sep. 2008; 
that claims priority to the following: Patent Cooperation 
Treaty Appl. No. PCT/US2006/036156, filed 18 Sep. 2006; 
U.S. Appl. No. 60/718, 393, filed 19 Sep. 2005; U.S. Appl. No. 
60/718,581, filed 19 Sep. 2005; U.S. Appl. No. 60/718,899, 
filed 20 Sep. 2005; U.S. Appl. No. 60/728,928, filed 20 Oct. 
2005; and U.S. Appl. No. 60/728,929 filed, 20 Oct. 2005; and 
U.S. application Ser. No. 10/351,664 filed, 27 Jan. 2003 now 
U.S. Pat. No. 8,041,436 that claims priority to U.S. Appl. No. 
60/373,404, filed 18 Apr. 2002. This application is also a 
continuation of U.S. application Ser. No. 12/1 95,353, filed on 
20 Aug. 2008, and issued as U.S. Pat. No. 8,180,464 on 15 
May 2012. 

This work was supported at least in part by NASA under 
NASA contract number GT3-52387. Accordingly, the United 
States government may have certain rights herein. 

TECHNICAL FIELD 

The subject matter described herein relates to controllers, 
systems, and methods for feedback control of various sys- 
tems. 


cations. The following brief synopsis is intended only to 
provide background on some exemplary applications to assist 
a person of ordinary skill in the art in understanding this 
disclosure more Hilly. Additional background information on 
5 a select number of exemplary industrial applications are also 
provided. 

BACKGROUND OF THE INVENTION 

io Classic Control Theory provides a number of techniques 
for an engineer to use in controller design. Existing control- 
lers for linear, time invariant, and single-input single-output 
plants are categorized in one taxonomy into three forms: the 
proportional/integral/derivative (PID) controllers, transfer 
15 function based (TFB) controllers, and state feedback (SF) 
controllers. The PID controller is defined by the equation: 

u=K p e+Kje+K D e Equation (1) 

where (u) is the control signal and (e) is the error between the 
20 set point and the process output (y) being controlled. As 
commonly used, the notation of a dot (•) above a variable 
indicates the use of a derivative of the variable where the order 
corresponds to the number of dots. 

The PID-type of controller has been employed in engineer- 
25 ing and other applications since the early 1920s. It is an 
error-based controller that does not require an explicit math- 
ematical model of the plant 110. 

Transfer Function Based Controller 

The transfer function based (TFB) controller is given in the 
30 form of 

U(s) = G c (s)E(s) Equation (2) 


BACKGROUND 

35 

A feedback (closed-loop) control system 10, as shown in 
Prior Art FIG. 1, is widely used to modify the behavior of a 
physical process, denoted as the plant 110, so it behaves in a 
specific desirable way over time. For example, it may be 
desirable to maintain the speed of a car on a highway as close 40 
as possible to 60 miles per hour in spite of possible hills or 
adverse wind; or it may be desirable to have an aircraft follow 
a desired altitude, heading and velocity profile independently 
of wind gusts; or it may be desirable to have the temperature 
and pressure in a reactor vessel in a chemical process plant 45 
maintained at desired levels. All of these industrial tasks are 
accomplished today by using traditional feedback control, 
and the above are examples of what automatic control sys- 
tems are designed to do, without human intervention. 

The key component in a feedback control system is the 50 
controller 120, which determines the difference between the 
output “y” of the plant 110, (e.g., the temperature or position) 
and its desired value, and produces a corresponding control 
signal “u” (e.g., modulating the power input to a heater or 
motor). 55 

The goal of controller design is to make the difference 
between the actual output of the plant and the desired output 
as small as possible as rapidly as possible. Today, controllers 
are employed in a large number of industrial control applica- 
tions and in areas like robotics, aeronautics, astronautics, 60 
motors, motion control and thermal control, just to name a 
few. 

INDUSTRIAL APPLICABILITY 


n(s ) 

gm = 

d[s) 


where U(s) and E(s) are Laplace Transforms of u and e 
defined above, and n(s) and d(s) are polynomials in s. The 
TFB controller can be designed using methods in control 
theory based on the transfer function model of the plant, 
Gp(s). A PID controller can be considered a special case of a 
TFB controller because it has an equivalent transfer function 
of s. 


k t 

G c (s) = k p + — +k d s 


Equation (3) 


State Feedback (SF) Controller 

The State Feedback (SF) controller can be defined by 

u=r+Kx Equation (4) 

where u is the control input, r is the setpoint for the output to 
follow and x is the state vector associated with the system, 
based on the state space model of the plant 110: 

x(t)=Ax(t)+Bu(t ) Equation (5) 

y(t) Cx(t)+Duil ) Equation (6) 

where x(t) is the state vector, y(t) is the output vector, u(t) is 
the input or control vector, A is the state matrix, B is the input 
matrix, C is the output matrix and D is the feed-through 
matrix. When the state x is not accessible, a state observer 
(SO): 


Hie present system and method has potential applicability 
to a wide range of different industrial and co mm ercial appli- 


x=Ax+Bu+L(y-y) 


Equation (7) 
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is usable to estimate the state x, where the (') symbol is used 
to denote an estimate or observed value of a given variable 
and L denotes the observer coefficient matrix. 

State Observers 

Observers extract real-time information of a plant’s 110 5 
internal state from its input-output data. The observer usually 
presumes precise model information of the plant 110, since 
performance is largely based on its mathematical accuracy. 
Most prior art closed loop controllers require both types of 
information. Such presumptions, however, often make the to 
method impractical in engineering applications, since the 
challenge for industry remains in constructing these models 
as part of the design process. Another level of complexity is 
added when gain scheduling and adaptive techniques are used 
to deal with nonlinearity and time variance, respectively. 15 
Disturbance Observers and Disturbance Rejection 

Recently, disturbance rejection techniques have been used 
to account for uncertainties in the real world and successfully 
control complex nonlinear systems. The premise is to solve 
the problem of model accuracy in reverse by modeling a 20 
system with an equivalent input disturbance d that represents 
any difference between the actual plant P and a derived/ 
selected model P„ of the plant, including external distur- 
bances w. An observer is then designed to estimate the dis- 
turbance in real time and provide feedback to cancel it. As a 25 
result, the augmented system acts like the model P„ at low 
frequencies, making the system behave like P„ and allowing a 
controller to be designed for P„. 

The most common of these techniques is the disturbance 
observer (DOB) structure (Endo, S., H. Kobayashi, C J. 30 
Kempf, S. Kobayashi, M. Tomizuka and Y. Hori (1996). 
“Robust Digital Tracking Controller Design for High-Speed 
Positioning Systems.” Control Eng. Practice. 4:4, 527-536; 
Kim, B. K„ H.-T. Choi, W. K. Chung and I. H. Suh (2002). 
“Analysis and Design of Robust Motion Controllers in the 35 
Unified Framework.” J. of Dynamic Systems, Measurement, 
and Control, 124, 313-321; Lee, H. S. and M. Tomizuka 
(1996). “Robust Motion Controller Design for High-Accu- 
racy Positioning Systems.” IEEE Trans. Ind. Electron.. 43:1, 
48-55; Tesfaye, A., H. S. Lee and M. Tomizuka (2000). “A 40 
Sensitivity Optimization Approach to Design of a Distur- 
bance Observer in Digital Motion Control.” IEEE/ASME 
Trans, on Mechatronics, 5:1, 32-38; Umeno, T. and Y. Hori 
(1991). “Robot Speed Control of DC Servomotors Using 
Modern Two Degrees of-Freedom Controller Design”. IEEE 45 
Trans. Ind. Electron., 38:5, 363-368). It uses a disturbance 
observer consisting of a binomial Q-filters and the inverse 
model of the plant. And the controller is called parameterized 
in the sense that the Q-filter is treated as a parameter. A model 
deliberately different from P is suggested in E. Schrijver and 50 
J. Van Dijk, “Disturbance Observers for Rigid Mechanical 
Systems: Equivalence, Stability, and Design,” Journal of 
Dynamic Systems, Measurement, and Control , vol. 1 24, no. 4, 
pp. 539-548, 2002 to facilitate design, but no guidelines are 
given other than it should be as simple as possible, while 55 
cautioning stability and performance may be in danger if the 
model is significantly different from the plant. Another 
obstacle is that a separate state observer must be designed to 
provide state feedback to the controller. In existing research, 
derivative approximates are used in this way but their effect 60 
on performance and stability is questionable. Furthermore, 
the controller design is dependent on the DOB design, mean- 
ing that derivative approximates can not be arbitrarily 
selected. 

Multiple DOBs were used to control a multivariable robot 65 
by treating it as a set of decoupled single-input single-output 
(SISO) systems, each with disturbances that included the 


4 

coupled dynamics (Bickel, R. andM. Tomizuka (1999). “Pas- 
sivity-Based Versus Disturbance Observer Based Robot Con- 
trol: Equivalence and Stability.”./ of Dynamic Systems, Mea- 
surement, and Control, 121, 41-47; Hori, Y., K. Shimura and 
M. Tomizuka (1992). “Position/Force Control of Multi-Axis 
Robot Manipulator Based on the TDOF Robust Servo Con- 
troller For Each Joint.” Proc. of ACC, 753-757; Kwon, S. J. 
and W. K. Chung (2002). “Robust Performance of the Mul- 
tiloop Perturbation Compensator.” IEEE/ASME Trans. 
Mechatronics, 7:2, 190-200; Schrijver, E. and J. Van Dijk 
(2002) Disturbance Observers for Rigid Mechanical Sys- 
tems: Equivalence, Stability, and Design.”./ of Dymamic Sys- 
tems, Measurement, and Control, 124, 539-548.) 

Another technique, referred to as the unknown input 
observer (UI0), estimates the states of both the plant and the 
disturbance by augmenting a linear plant model with a linear 
disturbance model (Burl, J. B. (1999). Linear Optimal Con- 
trol , pp. 308-314. Addison Wesley Longman, Inc., California; 
Franklin, G. F., J. D. Powell andM. Workman ( 1 998). Digital 
Control of Dynamic Systems , Third Edition, Addison Wesley 
Longman, California; Johnson, CD. (1971). “Accommoda- 
tion of External Disturbances in Linear Regulator and Servo- 
mechanism Problems.” IEEE Trans. Automatic Control, 
AC-1 6:6, 635-644; Liu, C-S., and H. Peng (2002). “Inverse- 
Dynamics Based State and Disturbance Observer for Linear 
Time-Invariant Systems.” J. of Dymamic Systems, Measure- 
ment, and Control, 124, 375-381; Profeta, J. A. Ill, W. G. Vogt 
and M. H. Mickle (1 990). “Disturbance Estimation and Com- 
pensation in Linear Systems.” IEEE Trans. Aerospace and 
Electronic Systems, 26:2, 225-231; Schrijver, E. and J. Van 
Dijk (2002) “Disturbance Observers for Rigid Mechanical 
Systems: Equivalence, Stability, and Design.” J. of Dynamic 
Systems, Measurement, and Control, 124, 539-548). Unlike 
the DOB structure, the controller and observer can be 
designed independently, like a Luenberger observer. How- 
ever, it still relies on a good mathematical model and a design 
procedure to determine observer gains. An external distur- 
bance w is generally modeled using cascaded integrators 
(1/s"). When they are assumed to be piece-wise constant, the 
observer is simply extended by one state and still demon- 
strates a high degree of performance. 

Extended State Observer (ESO) 

In this regard, the extended state observer (ESO) is quite 
different. Originally proposed by Han, J. (1999). “Nonlinear 
Design Methods for Control Systems.” Proc. lAth IFAC 
World Congress, in the form of a nonlinear UIO and later 
simplified to a linear version with one tuning parameter by 
Gao, Z. (2003). “Scaling and Parameterization Based Con- 
troller Tuning.” Proc. of ACC, 4989-4996, the ESO combines 
the state and disturbance estimation power of a UIO with the 
simplicity of a DOB. One finds a decisive shift in the under- 
lying design concept as well . The traditional observer is based 
on a linear time-invariant model that often describes a non- 
linear time-varying process. Although the DOB and UIO 
reject input disturbances for such nominal plants, they leave 
the question of dynamic uncertainty mostly unanswered. The 
ESO, on the other hand, addresses both issues in one simple 
framework by formulating the simplest possible design 
model P rf =l/s" for a large class of uncertain systems. The 
design model P rf is selected to simplify controller and 
observer design, forcing P to behave like the design model P rf 
at low frequencies rather than P„. As a result, the effects of 
most plant dynamics and external disturbances are concen- 
trated into a single unknown quantity. The ESO estimates this 
quantity along with derivatives of the output, giving way to 
the straightforward design of a high performance, robust, 
easy to use and affordable industrial controller. 
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Active Disturbance Rejection Control (ADRC) 

Originally proposed by Han, J. (1999). “Nonlinear Design 
Methods for Control Systems.” Proc. 14 th IFAC World Con- 
gress, a nonlinear, non-parameterized active disturbance 
rejection control (ADRC) is a method that uses an ESO. A 5 
linear version of the ADRC controller and ESO were param- 
eterized for transparent tuning by Gao, Z. (2003). “Scaling 
and Parameterization Based Controller Tuning.” Proc. of 
ACC, 4989-4996. Its practical usefulness is seen in a number 
of benchmark applications already implemented throughout to 
industry with promising results (Gao, Z., S. Hu and F. Jiang 
(2001). “A Novel Motion Control Design Approach Based on 
Active Disturbance Rejection.” Proc. of 40th IEEE Confer- 
ence on Decision and Control ; Goforth, F. (2004). “On 
Motion Control Design and Timing Techniques.” Proc. of 15 
ACC, Hu, S. (2001). “On High Performance Servo Control 
Solutions for Hard Disk Drive.” Doctoral Dissertation, 
Department of Electrical and Computer Engineering, Cleve- 
land State University; Hou, Y., Z. Gao, F. Jiang and B. T. 
Boulter (2001). “Active Disturbance Rejection Control for 20 
Web Tension Regulation.” Proc. of 40 th IEEE Conf. on Deci- 
sion and Control, Huang, Y., K. Xu and J. Han (2001 ). “Flight 
Control Design Using Extended State Observer and Nons- 
mooth Feedback.” Proc. of 40 th IEEE Conf. on Decision and 
Control, Sun, B and Z. Gao (2004). “A DSP-Based Active 25 
Disturbance Rejection Control Design for a 1 KW H-Bridge 
DC-DC Power Converter.” IEEE Trans, on Ind. Electronics-, 
Xia, Y., L. Wu, K. Xu, and J. Han (2004). “Active Disturbance 
Rejection Control for Uncertain Multivariable Systems With 
Time-Delay, 2004 Chinese Control Conference. It also was 30 
applied to a fairly complex multivariable aircraft control 
problem (Huang, Y., K. Xu and J. Han (2001 ). “Flight Control 
Design Using Extended State Observer and Nonsmooth 
Feedback. Proc. of 40 th IEEE Conf. on Decision and Con- 
trol). 35 

What is needed is a control framework for application to 
systems throughout industry that are complex and largely 
unknown to the personnel often responsible for controlling 
them. In the absence of required expertise, less tuning param- 
eters are needed than current approaches, such as multi-loop 40 
PID, while maintaining or even improving performance and 
robustness. 

INDUSTRIAL APPLICATION 

45 

A selection of industrial applications is presented herein to 
provide a background on some selected applications for the 
present extended active disturbance rejection controller. 
Tracking Control 

Tracking Control refers to the output of a controlled system 50 
meeting design requirements when a specified reference tra- 
jectory is applied. Oftentimes, it refers to how closely the 
output y compares to the reference input r at any given point 
in time. This measurement is known as the error e=r-y. 

Control problems are categorized in two major groups; 55 
point-to-point control and tracking control. Point-to-point 
applications usually call for a smooth step response with 
small overshoot and zero steady state error, such as when 
controlling linear motion from one position to the next and 
then stopping. Since the importance is placed on destination 60 
accuracy and not on the trajectory between points, conven- 
tional design methods produce a bandwidth limited controller 
with inherent phase lag in order to produce a smooth output. 
Tracking applications require a good tracking of a reference 
input by keeping the error as small as possible at all times, not 65 
just at the destination. This is particularly co mm on in con- 
trolling a process that never stops. Since the importance is 
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placed on accurately following a changing reference trajec- 
tory between points, phase lags are avoided as much as pos- 
sible because they may lead to significant errors in the tran- 
sient response, which lasts for the duration of the process in 
many applications. A common method of reducing phase lag 
is to increase the bandwidth, at the cost of increased control 
effort (energy) and decreased stability margin. 

Various methods have been used to remove phase lag from 
conventional control systems. All of them essentially modify 
the control law to create a desired closed loop transfer func- 
tion equal to one. As a result, the output tracks the reference 
input without much phase lag and the effective bandwidth of 
the overall system is improved. The most common method is 
model inversion where the inverse of the desired closed loop 
transfer function is added as a prefilter. Another method pro- 
posed a Zero Phase Error Tracking Controller (ZPETC) that 
cancels poles and stable zeros of the closed loop system and 
compensates for phase error introduced by un-cancelable 
zeros. Although it is referred to as a tracking controller, it is 
really a prefilter that reduces to the inverse of the desired 
closed loop transfer function when unstable zeros are not 
present. Other methods consist of a single tracking control 
law with feed forward terms in place of the conventional 
feedback controller and prefilter, but they are application 
specific. However, all of these and other previous methods 
apply to systems where the model is known. 

Model inaccuracy can also create tracking problems. The 
performance of model-based controllers is largely dependent 
on the accuracy of the model. When linear time-invariant 
(LTI) models are used to characterize nonlinear time-varying 
(NTV) systems, the information becomes inaccurate over 
time. As a result, gain scheduling and adaptive techniques are 
developed to deal with nonlinearity and time variance, 
respectively. However, the complexity added to the design 
process leads to an impractical solution for industry because 
of the time and level of expertise involved in constructing 
accurate mathematical models and designing, tuning, and 
maintaining each control system. 

There have been a number of high perfomiance tracking 
algorithms that consist of three primary components: distur- 
bance rejection, feedback control, and phase lag compensa- 
tion implemented as a prefilter. First, disturbance rejection 
techniques are applied to eliminate model inaccuracy with an 
inner feedback loop. Next, a stabilizing controller is con- 
structed based on a nominal model and implemented in an 
outer feedback loop. Finally, the inverse of the desired closed 
loop transfer function is added as a prefilter to eliminate phase 
lag. Many studies have concentrated on unifying the distur- 
bance rejection and control part, but not on combining the 
control and phase lag compensation part, such as the RIC 
framework. Internal model control (IMC) cancels an equiva- 
lent output disturbance. B. Francis and W. Wonham, “The 
Internal Model Principal of Control Theory Automatica, vol 
12, 1976, pp. 457-465. E. Schrijver and J. Van Dijk, “Distur- 
bance Observers for Rigid Mechanical Systems: Equiva- 
lence, Stability, and Design,” Journal of Dynamic Systems, 
Measurement, and Control, vol. 124, December 2002, pp. 
539-548 uses a basic tracking controller with a DOB to con- 
trol a multivariable robot. The ZPETC has been widely used 
in combination with the DOB framework and model based 
controllers. 

Thus, having reviewed prior art in tracking control, the 
application now describes example systems and methods of 
tracking control employing predictive ADRC. 

Web Processing Applications 

Web tension regulation is a challenging industrial control 
problem. Many types of material, such as paper, plastic film, 
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cloth fabrics, and even strip steel are manufactured or pro- 
cessed in a web form. Hie quality of the end product is often 
greatly affected by the web tension, making it a crucial vari- 
able for feedback control design, together with the velocities 
at the various stages in the manufacturing process. The ever- 5 
increasing demands on the quality and efficiency in industry 
motivate researchers and engineers alike to explore better 
methods for tension and velocity control. However, the highly 
nonlinear nature of the web handling process and changes in 
operating conditions (temperature, humidity, machine wear, to 
and variations in raw materials) make the control problem 
challenging. 

Accumulators in web processing lines are important ele- 
ments in web handling machines as they are primarily respon- 
sible for continuous operation of web processing lines. For 15 
this reason, the study and control of accumulator dynamics is 
an important concern that involves a particular class of prob- 
lems. The characteristics of an accumulator and its operation 
as well as the dynamic behavior and control of the accumu- 
lator carriage, web spans, and tension are known in the art. 20 

Both open-loop and closed-loop methods are commonly 
used in web processing industries for tension control pur- 
poses. In the open-loop control case, the tension in a web span 
is controlled indirectly by regulating the velocities of the 
rollers at either end of the web span. An inherent drawback of 25 
this method is its dependency on an accurate mathematical 
model between the velocities and tension, which is highly 
nonlinear and highly sensitive to velocity variations. Never- 
theless, simplicity of the controller outweighs this drawback 
in many applications. Closing the tension loop with tension 30 
feedback is an obvious solution to improve accuracy and to 
reduce sensitivity to modeling errors. It requires tension mea- 
surement, for example, through a load cell, but is typically 
justified by the resulting improvements in tension regulation. 

Most control systems will unavoidably encounter distur- 35 
bances, both internal and external, and such disturbances 
have been the obstacles to the development of high perfor- 
mance controller. This is particularly true for tension control 
applications and, therefore, a good tension regulation scheme 
must be able to deal with unknown disturbances. In particular, 40 
tension dynamics are highly nonlinear and sensitive to veloc- 
ity variations. Further, process control variables are highly 
dependent on the operating conditions and web material char- 
acteristics. Thus, what are needed are systems and methods 
for control that are not only overly dependent on the accuracy 45 
of the plant model, but also suitable for the rejection of sig- 
nificant internal and external disturbances. 

Jet Engine Control Applications 

A great deal of research has been conducted towards the 
application of modern multivariable control techniques on 50 
aircraft engines. The majority of this research has been to 
control the engine at a single operating point. Among these 
methods are a multivariable integrator windup protection 
scheme (Watts, S. R. and S. Garg (1996). “An Optimized 
Integrator Windup Protection Technique Applied to a Turbo- 55 
fan Engine Control,” AIAA Guidance Navigation and Control 
Conf.), a tracking filter and a control mode selection for 
model based control (Adibhatla S. and Z. Gastineau (1994). 
“Tracking Filter Selection And Control Mode Selection For 
Model Based Control.” AIAA 30th Joint Propulsion Confer- 60 
ence and Exhibit), an H m method and linear quadratic Gaus- 
sian with loop transfer recovery method (Watts, S. R. and S. 
Garg (1995).” A Comparison Of Multivariable Control 
Design Techniques For A Turbofan Engine Control.” Inter- 
national Gas Turbine and Aeroengine Congress and Expo.), 65 
and a performance seeking control method (Adibhatla, S. and 
K. L. Johnson (1993). “Evaluation of a Nonlinear Psc Algo- 
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rithm on a Variable Cycle Engine.” AIAA/SAE/ASME/ASEE 
29th Joint Propulsion Conference and Exhibit). Various 
schemes have been developed to reduce gain scheduling 
(Garg, S. (1997). “A Simplified Scheme for Scheduling Mul- 
tivariable Controllers.” IEEE Control Systems) and have even 
been combined with integrator windup protection and H 2 
(Frederick, D. K., S. Garg and S. Adibhatla (2000). “Turbofan 
Engine Control Design Using Robust Multivariable Control 
Technologies. IEEE Trans, on Control Systems Technology >). 

Conventionally, there have been a limited number of con- 
trol techniques for full flight operation (Garg, S. (1997). “A 
Simplified Scheme for Scheduling Multivariable Control- 
lers.” IEEE Control Systems', and Polley, J. A., S. Adibhatla 
and P J. Hoffman (1988). “Multivariable Turbofan Engine 
Control for Full Conference on Decision and Control Flight 
Operation.” Gas Turbine and Expo). However, there has been 
no development of tuning a controller for satisfactory perfor- 
mance when applied to an engine. Generally, at any given 
operating point, models can become inaccurate from one 
engine to another. This inaccuracy increases with model com- 
plexity, and subsequently design and tuning complexity. As a 
result, very few of these or similar aircraft design studies have 
led to implementation on an operational vehicle. 

The current method for controlling high performance jet 
engines remains multivariable proportional -integral (PI) con- 
trol (Edmunds, J. M. (1979). “Control System Design Using 
Closed-Loop Nyquist and Bode Arrays.” Int. J. on Control. 
30:5, 773-802, and Polley, J. A., S. Adibhatla and P J. Hoff- 
man (1 988). “Multivariable Turbofan Engine Control for Full 
Conference on Decision and Control. Flight Operation.” Gas 
Turbine and Expo). Although the controller is designed by 
implementing Bode and Nyquist techniques and is tunable, a 
problem remains due to the sheer number of tuning param- 
eters compounded by scheduling. 

SUMMARY OF INVENTION 

The present system, method, process, and device for con- 
trolling systems is applicable for a variety of different appli- 
cation and has utility in a number of different industrial fields. 
The following brief summary section is provided to facilitate 
a basic understanding of the nature and capabilities of the 
present system and method. This summary section is not an 
extensive nor comprehensive overview and is not intended to 
identify key or critical elements of the present systems or 
methods or to delineate the scope of these items. Rather this 
brief summary is intended to provide a conceptual introduc- 
tion in simplified form as an introduction for the more 
detailed description presented later in this document. 

The present application describes selected embodiments of 
the ADRC controller that comprise one or more computer 
components, that extend, build upon and enhance the use of 
ADRC controllers and provide enhanced performance and 
utility. Specifically, in one aspect the ADRC controller uti- 
lizes a predictive computer component or module that in one 
embodiment predicts future values of the plant output and in 
a second embodiment predicts future estimates of the system 
state and generalized disturbance. The generalized distur- 
bance includes dynamics of the plant itself such that the plant 
is effectively reduced to a cascaded integral plant. In the 
various embodiments a variety of estimates are available and 
introduced. 

In still another aspect, the ADRC controller further 
includes a model of the plant dynamics. The model is used in 
one embodiment to improve the state estimator or predictor. 
In a second embodiment the model is used to provide an 
enhanced control law that provides improved performance. In 
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both aspects of the embodiments discussed in this paragraph, 
the generalized disturbance captures model errors and dis- 
crepancies and allows the system to eliminate errors caused 
by model errors, discrepancies, and assumptions. In the mul- 
tiple embodiments, a variety of models are available for use 
by one of ordinary skill in the art to achieve the desired 
performance relative to the type of plant being controlled. 

In yet another aspect, the ADRC controller further includes 
a non-linear, or discrete time optimal control law. The non- 
linear control law of the present embodiment utilizes non- 
linear control laws to improve the overall performance of an 
ADRC controller as compared to an ADRC controller with a 
linear proportional -derivative based control law. 

In still another aspect, additional knowledge of derivative 
plant output gathered from high quality sensor information or 
direct measurement sensors (e.g. velocity and acceleration 
sensors) is used to reduce the order of the system state and 
disturbance estimator. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The accompanying figures depict multiple embodiments of 
a device and method for the separation of constituents from a 
flow. A brief description of each figure is provided below. 
Elements with the same reference numbers in each figure 
indicate identical or functionally similar elements. Addition- 
ally, the left-most digit(s) of a reference number identifies the 
drawings in which the reference number first appears. 

FIG. 1 is a block diagram of a Prior Art feedback controller. 

FIG. 2 is a block diagram of an ADRC configuration for a 
2nd order plant. 

FIG. 3 is a block diagram of a predictive ADRC configu- 
ration. 

FIG. 4 is a block diagram of a 2 degree of freedom transfer 
function based (TFB) control structure. 

FIG. 5 is a block diagram of an implementation of a 
PADRC configuration based on output prediction. 

The present invention is described with reference to block 
diagrams and operational flow charts. It is to be understood 
that the functions/acts noted in the blocks may occur out of 
the order noted in the operational illustrations. For example, 
two blocks shown in succession may in fact be executed 
substantially concurrently or the blocks may sometimes be 
executed in the reverse order, depending upon the function- 
ality/acts involved, including for example executing as asyn- 
chronous threads on a processor. Although some of the dia- 
grams include arrows on communication paths to show a 
primary direction of communication, it is to be understood 
that communication may occur in the opposite direction to the 
depicted arrows. 

DETAILED DESCRIPTION 

Multiple embodiments of a system, device, and method for 
the control of systems are presented herein. Those of ordinary 
skill in the art can readily use this disclosure to create alter- 
native embodiments using the teaching contained herein. 
Lexicon 

The following terms used herein have the meanings as 
follows. 

As used herein the term “computer component” refers to a 
computer and elements of a computer, such as hardware, 
firmware, software, a combination thereof, or software in 
execution. For example, a computer component can include 
by way of example, a process running on a processor, a 
processor, an object, an executable, an execution thread, a 
program, and a computer itself. One or more computer com- 
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ponents can in various embodiments reside on a server and 
the server can be comprised of multiple computer compo- 
nents. One or more computer components are in some cases 
referred to as computer systems whereby one or more com- 
5 puter components operate together to achieve some function- 
ality. One or more computer components can reside within a 
process and/or thread of execution and a computer compo- 
nent can be localized on one computer and/or distributed 
between two or more computers, 
to “Software”, as used herein, includes but is not limited to, 
one or more computer readable and/or executable instruc- 
tions that cause a computer or other electronic device to 
perform functions, actions and/or behave in a desired maimer. 
The instructions may be embodied in various forms like rou- 
ts tines, algorithms, modules, methods, threads, and/or pro- 
grams. Software may also be implemented in a variety of 
executable and/or loadable forms including, but not limited 
to, a stand-alone program, a function call (local and/or 
remote), a servelet, an applet, instructions stored in a memory, 
20 part of an operating system or browser, and the like. It is to be 
appreciated that the computer readable and/or executable 
instructions can be located in one computer component and/ 
or distributed between two or more communicating, co-op- 
erating, and/or parallel processing computer components and 
25 thus can be loaded and/or executed in serial, parallel, mas- 
sively parallel and other maimers. It will be appreciated by 
one of ordinary skill in the art that the form of software may 
be dependent on, for example, requirements of a desired 
application, the environment in which it runs, and/or the 
30 desires of a designer/programmer or the like. 

“Computer communications”, as used herein, refers to a 
communication between two or more computers and can be, 
for example, a network transfer, a file transfer, an applet 
transfer, an email, a hypertext transfer protocol (HTTP) mes- 
35 sage, a datagram, an object transfer, a binary large object 
(BLOB) transfer, and so on. A computer communication can 
occur across, for example, a wireless system (e.g., IEEE 
802.1 1 ), an Ethernet system (e.g., IEEE 802.3), a token ring 
system (e.g., IEEE 802.5), a local area network (LAN), a wide 
40 area network (WAN), a point-to-point system, a circuit 
switching system, a packet switching system, and so on. 

.An “operable connection” is one in which signals and/or 
actual communication flow and/or logical communication 
flow may be sent and/or received. Usually, an operable con- 
45 nection includes a physical interface, an electrical interface, 
and/or a data interface, but it is to be noted that an operable 
connection may consist of differing combinations of these or 
other types of comiections sufficient to allow operable con- 
trol. 

50 As used herein, the term “signal” may take the form of a 
continuous waveform and/or discrete value(s), such as digital 
value(s) in a memory or register, present in electrical, optical 
or other form. 

The term “controller” as used herein indicates a method, 
55 process, or computer component adapted to control a plant 
(i.e. the system to be controlled) to achieve certain desired 
goals and objectives. 

As used herein the following symbols are used to describe 
specific characteristics of the variable. A variable with a 
60 overtop the variable label, unless other indicated to the con- 
trary, indicates that the variable is an estimate of the actual 
variable of the same label. A variable with a “” overtop the 
variable label, unless otherwise indicated, indicates that the 
variable is an nth derivative of that variable, where n equals 
65 the number of dots. 

To the extent that the term “includes” is employed in the 
detailed description or the claims, it is intended to be inclusive 
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in a manner similar to the term “comprising” as that term is 
interpreted when employed as a transitional word in a claim. 

To the extent that the term “or” is employed in the claims 
(e.g., A or B) it is intended to mean “A or B or both”. When the 
author intends to indicate “only A or B but not both”, then the 
author will employ the term “A or B but not both”. Thus, use 
of the term “or” in the claims is the inclusive, and not the 
exclusive, use. 

Active Disturbance Rejection Control (ADRC) 

Active Disturbance Rejection Control (“ADRC”) is a 
method or process and, when embedded within a computer 
system, a device (i.e., a controller that processes information 
and signals using the ADRC algorithms), that effectively 
rejects both unknown internal dynamics of a system of an 
arbitrary order to be controlled, or plant 110, and external 
disturbances in real time without requiring detailed knowl- 
edge of the plant 110 dynamics, which is required by most 
existing control design methods. 

The only information needed to configure an ADRC con- 
troller is knowledge of the relative order of the plant 110 and 
the high frequency gain of the plant 110. In comparison, 
traditional model-based control methods require accurate 
mathematical models to conform with prevailing model- 
based design methods. In some embodiments, the ADRC 
controller leverages knowledge gained from model-based 
controller designs to improve control performance. The 
ADRC approach is unique in that it first forces an otherwise 
nonlinear, time- varying and uncertain plant 110 to behave as 
a simple linear cascade integral plant that can be easily con- 
trolled. To this end, the unknown, nonlinear, and time- varying 
internal dynamics is combined with the external disturbances 
to form what is denoted as the generalized disturbance, which 
is then estimated and rejected, thus greatly simplifying the 
control design. 

The embodiments of the ADRC controller design are 
founded on the basic theory that uncertainties in dynamics 
and unknown external disturbances in the operating environ- 
ment are a key challenge to be dealt with in control engineer- 
ing practice. Traditional control systems are incapable of 
providing a solution to controlling and mitigating the impact 
of these uncertainties. The ADRC controller in contrast is 
specifically adapted to mitigate and eliminate the impact of 
these uncertainties, thus providing a dynamically simplified 
system for control. 

One previous major characteristic of the ADRC controller 
is the parameterization of the controller as a function of 
control loop bandwidth. This parameterization makes all of 
the ADRC controller parameters functions of the control loop 
bandwidth; see, U.S. patent application Ser. No. 10/351,664, 
titled, “Scaling and Parameterizing A C ontroller” to Zhiqiang 
Gao, which is hereby incorporated by reference in its entirety. 

An embodiment of a parameterized ADRC controller 
applied to a second-order plant 202 is illustrated in FIG. 2. 
The system to be controlled is a second-order plant 202 and 
represents an arbitrary plant 110 with predominantly second- 
order dynamics. The second-order plant 202 has an output (y) 
in response to an input that comprises a control signal (u) with 
superimposed external disturbances (d). 

Consider a general uncertain, nonlinear and time-varying 
second-order plant 202 of the form: 

y~j{y,y, d,t)+bu Equation (8) 

where y is system output, u is the control signal, b is a 
constant, d is external input disturbance, f(y,y,d,t) is treated as 
the generalized disturbance. When this generalized distur- 
bance is cancelled, the system is reduced to a simple double- 
integral plant with a scaling factor b, which can be easily 
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controlled. Here we assume that the approximate value of b, 
denoted as b, is given or estimated. As shown in FIG. 2, the 
scaling factor (b) is used to scale 208, the control output (uft). 
In other embodiments, the scaling factor is implicit within the 
5 linear PD control law 206. 

To capture the information of the generalized disturbance 
f(y,y,d,t) and cancel it from the system dynamics leads us to 
the Extended State Observer (ESO) 204, which is now pre- 
sented. Let x 1 =y, x 2 =y, and x 3 =f, the description of the above 
to plant 202 is rewritten in state space form with f(y,y,d,t), or 
simply f, treated as an additional state: 

(x = A x x + B x u + Et Equation (9) 

15 i, = G, 

where 
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The corresponding state observer is 

I i = A x z + B z u + Uy m - C X Z ) Equation (10) 

30 | f = Vz 

where y„, is the measured system output, B_=[0 b 0] r , V=[0 0 
1J. and L=[lj 1, 1 3 ] 7 is the observer gain vector. With the 
35 parameterization technique as described more fully in U.S. 
patent application Ser. No. 10/351,664, L is obtained by solv- 
ing the equation Z(s)=lsI-A 3 .+LC x l=(s+co c) ) 3 ; that is, all the 
eigenvalues of the observer are placed at -tu 0 with 

/i=3oj, j ./ 2 = 3oj,/// 3 =oj 0 3 Equation (11) 

40 

With a well-tuned extended state observer 204, z 15 z 2 , and 
z 3 will closely track y, y, and f respectively, i.e. z^y, z 2 «y, and 
z 3 ~f 

Using the information provided by the extended state 
4 5 observer 204, the generalized disturbance f is rejected by the 
control law 

u=(u 0 -z 3 )/b Equation ( 1 2) 

which reduces the plant 110 to a unity gain, double integral 

50 P lant 

y=u 0 Equation (13) 

that can be easily controlled with linear proportional -deriva- 
tive (PD) with feedforward 

55 u 0 =g p J,e,e,(ti c )+r Equation (14) 

where g pd (e,e,a) c ) is the parameterized linear PD controller 
with the loop bandwidth of w c 

gpjte.e,™ J=(<o/e+2(o c e Equation (15) 

60 with the error term defined as 

e=r-z l e =r-z 2 Equation (16) 

where the desired trajectory (r) is provided by the user. 

Note that the linear PD control law 206 in Equation (14) 
65 corresponds to the co mm on PD gains of 

k r < 11 / and 2oq 


Equation (17) 
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Using these details, the ADRC controller is generalized in 
its various embodiments to control an arbitrary n rf '-order sys- 
tem y ( " ) =f(y,y, . . . , y ( ” _1) ,d,t)+bu, in which the external 
disturbances and unknown internal dynamics are rejected. In 
a similar manner, the plant 110 is reduced to a cascaded 5 
integral plant y (,,) =u, which is then controlled by a general- 
ized proportional -derivative (PD) control law, also described 
as a linear PD control law 206 similar to Equation (14). 

Use of Enhanced Sensor Information to Improve Control 10 
Performance 

In other embodiments of the ADRC controller, the linear 
PD control law 206 is updated to enhance performance by 
using direct measurement of the system output states (y m ) 
instead of the estimated states (z t , z 2 ). Specifically, when the 15 
measured output y m (t) is sufficiently noise free and y m (t) is 
obtainable, either via direct measurement or from y m (t), the 
control law in Equation (14) can be implemented with 

e=r-y m e=r-y m Equation (18) ^ 

That is, the output measurement and its derivative are used 
in the control law, instead of their estimates, z l and z 2 . Simi- 
larly for the n*-order plant 110, the output measurement y„, 
and its derivatives from the Uto the (n- 1 ) rf '-order available, 
can also be used directly in the PD control law, also referred 
to as the linear PD control law 206, instead of their estimates 
from the Extended State Observer (ESO) 204. Using the 
direct measurement of the output measurement state y m and 30 
its derivatives reduces the phase lag and improve control 
system performance, if the output measurement and its 
derivatives are not corrupted by noise. 

Multiple Extended States 35 

The Extended State Observer 204 in alternative embodi- 
ments includes more than one extended states, such as 
embodiments where the 2nd order plant 202 is replaced with 
a higher order plant 110. In one embodiment, when aug- 40 
mented with two extended states, x 3 =f and x 4 =f, the extended 
plant state space model can be written as follows. 

tx = A x x + B x u + E? Equation (19) 45 

\y = C x x 


where 
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With z u z 2 , z 3 , and z 4 tracking x l5 x 2 , x 3 , and x 4 respec- 
tively, the 4 2, '-order ESO 204 is employed to estimate f. 


( Z = A X Z + B z u + Uy m - c x z ) Equation (20) 

\f = Vz 


where B, |0 b 0 Ofand V=[0 0 10], 

In this manner, any embodiment of an ADRC controller is 
adaptable for use with an arbitrary number (h) extended state 
observer states with is generalized to control an n f, ’-order 
system. 

Discrete Implementation 

ADRC algorithms are implemented in discrete time, with 
ESO 204 discretized in both Euler and ZOE1, and in both the 
current estimator and predictive estimator form. 

In discrete implementation, the ESO 204 pole location is 
determined in the Z domain by solving the equation X(z)=l zl- 
0 +L p C=(z-flf, where |3=e-“-’ :r . 

Predictive Discrete Extended State Observer (PDESO) is 
in predictive estimator form, as shown in the following equa- 
tion. 


I z[k + 1] = 4>z[A] + E u[k] + L p (y m [k ] - zi [A]) Equation (21) 

U[A] = Z[A] 

Current Discrete Extended State Observer (CDESO) is in 
current estimator form, as shown in the following equation. 

tz [ A + 1] = 4>z[A] + r«[A] + Lp(y m [k\ - zi [A]) Equation (22) 

lz[A] = z[A] + UU m [A]-Zi[A]) 

The gain matrices of different ADRC discrete implemen- 
tations are shown in the following Table I. 

TABLE I 


PDESO AND CDESO COEFFICIENT MATRICES 


Type <E> T L p L c 


Euler 

1 T O' 


0 ' 


3(1 -/?) 


l-/? 3 


0 1 T 


bT 


3(1 -/?) 2 /T 


(l-/?) 2 (2 + / S)/T 


. 001 . 


. 0 . 


ji-/J) 3 /t 2 


(1 -/?) 3 /t 2 


T 2 

1 T T 

0 1 T 
.0 0 1 


r.T 2 1 


3(1 -P) 


1 -0 1 


b 2 
bT 
0 


3(l-/?) 2 (5+/?)/2T 

(1-/J) 3 /T 2 


3(1 -/?) 2 (1 +/?)/2T 
(l-/?) 3 /T 2 
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This discrete implementation of ADRC is also generaliz- 
able to n th -order plants 110. 

Predictive Active Disturbance Rejection Controller 

The thrust of the many of the embodiments presented here 
is the addition of the predictive functions to extend the capa- 
bilities of the original components of the ADRC. This extends 
the ADRC controller and generates embodiments of predic- 
tive versions of the ADRC, such as the predictive ADRC 300, 
depicted in block diagram form in FIG. 3 In this embodiment, 
the linear PD control law 206 in the predictive ADRC 300 is 
now a predictive control law 306. Similarly, the observer, or 
more specifically the extended state observer 204 from the 
traditional ADRC 200 is now embodied as a predictive state 
and disturbance observer 304. In one embodiment the predic- 
tive state and disturbance observer 304 is a predictive 
extended state observer (PESO) where the prediction occurs 
within the bounds of the ESO 204 to where the estimated 
states z„(t) are predicted to form the predictive state and 
disturbance observer 304, and outputs the prediction of the 
estimated states z„(t+x). In a second embodiment, a state 
predictor method is used to extend the design of the ESO 204 
to estimate the states z„(t) and form the predictive state and 
disturbance observer 304. 

In other embodiments, other estimators and predictors, 
either in an integrated form or in combination are used to 
generate similar information as that generated by the predic- 
tive state and disturbance observer 304 and used by the 
embodiments predictive ADRC 300. One of these other 
embodiments adopts a disturbance observer (DOB) estimator 
structure to estimate the disturbance and a state observer and 
a predictor element that together form the predictive state and 
disturbance observer 304. In a second other embodiment, the 
predictive state and disturbance observer 304 utilizes an 
unknown input observer (UIO) estimator. Other estimators, 
including combinations of estimators such as the combina- 
tion of a state observer and a DOB, are suitable for estimation 
of both the plant 110 output and the known part of the general 
disturbance (f„) are adaptable using the disclosure provided 
herein to those of ordinary skill in the art to serve as a pre- 
dictive state and disturbance observer 304. 

The challenge in constructing the various embodiments of 
a predictive controller within the general purpose ADRC 
controller architecture is how the various signals used by the 
ADRC algorithm are predicted and used. The signals used by 
the ADRC algorithm originate from a multitude of sources in 
the various embodiments, including in one embodiment the 
extended state observer 204 (ESO). In another embodiment, a 
simplified, reduced-order ESO (RESO) is used when either: 
(i) quality measurements of the output signals allows numeri- 
cal differentiation to estimate derivatives of the output signal; 
or (ii) direct measurement of derivatives are possible (e.g. 
velocity and acceleration sensors), are available. 

In another embodiment, a nonlinear proportional-deriva- 
tive (NPD) control law is used to further optimize the perfor- 
mance of the present ADRC system. The RESO aspect of the 
ADRC controller provides alternative embodiments of ESO 
204 and RESO; a general form of RESO with multiple 
extended states for an lU-order plant 110, and incorporation 
of the known plant 110 dynamics into ESO 204 and control 
law. 

Furthermore, to facilitate practical implementations of the 
new techniques, we have provided at selected locations 
herein, the transfer function equivalent predictive ADRC 400, 
with H(s) 402 and G^s) 404, in a two-degree-of-freedom 
(2dof) structure shown in FIG. 4. Similarly, for a transfer 
function equivalent predictive ADRC 400, the plant dynamics 
are represented as G^(s) 402. For various forms of the Pre- 
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dictive ADRC 400 described herein, these transfer function 
expressions, together with the NPD and RESO, are novel and 
significant even when the amount of prediction is set to zero. 
All of these new developments help to make the predictive 
5 ADRC 300 possible and to provide the means of implemen- 
tation the present embodiments in a diverse range of indus- 
trial machinery, and they are illustrated in details in the fol- 
lowing sections. 

Basic Design Principles for Predictive ADRC 
to In the real world, the sensor feedback, the generation of the 
control signal and its arrival at the actuator, and the estimation 
of the states and disturbances all contain time delays and 
phase lags, which can endanger system stability. To accom- 
modate such delays and lags, we propose the following Pre- 
15 dictive ADRC control law: 

u(f) = (g(.e(t+! e ),e(t+l i ), < a c )+r(t+l u )-Xt+lJ))/b Equation (23) 

where g(e,e,co £7 ) is the parameterized linear PD controller 206, 
or nonlinear proportional derivative (NPD) controller, Equa- 
20 tion (68), l e , 4, 1„, and l^are positive real numbers representing 
the prediction horizons, a special case of which is 
l e =4=l„=4=T and the corresponding control law is 

u(t)=(g(e(t+x),e(t+x),u> c )+f(t+x)-Xt+xyb Equation (24) 

Also, if a part of the plant 110 dynamics, f\ is given, it can 
be incorporated into ADRC, as shown later in this disclosure. 
In the context of Predictive ADRC 300, the knowledge of f„ 
can be used to help predict its future value, f„(t+4). The 
. corresponding control law for Equation (23) and Equation 
(24) are 

uit)=igie(t+l t ,),e(t+l i ),ii^ c )+r(t+lJ-f n (t+lj- n )-fJt+lj- u ))/b Equation (25) 


u(t)=(g(e(t+x),e(t+l- z ),o> c )+Ht+x)-f n (t+x)-f u (t+x))/b Equation (26) 

respectively. 

Described below are various methods to obtain the predic- 
411 tion of e(t), e(t), and f(t), respectively, so that these predictive 
ADRC 300 control laws can be implemented. 

Prediction of e and e 

Two embodiments of methods for predicting the values of 
error e(t) and its derivative are presented herein. 

Method 1 : Using Taylor Series, the prediction of e(t), e(t) 
can be obtained as e(t+l e )=e(t)+e(t)l e or 


50 e(t + /„) = eU) + e(t)l, + 

and e(t+4)=e(t)+e(t)4. Here e(t) can be obtained directly 
from r-y, or from r-z 2 using regular ESO, or i:-z 1 using 
55 RESO. Here e(t) can be obtained as e(t)=f(t)-y(t)=r(t)-(z 3 
(t)+bu(t)) with ESO or e(t)=r(t)-y(t)“f(t)-(z 2 (t)+bu(t)) with 
RESO. 

Method 2: Since e(t+x)=r(t+x)-z 1 (t+x) and e(t+x)=r(t+x)- 
z 2 (t+x), e(t+x) and e(t+x) can be obtained directly using the 
60 fiiture reference signals r(t+x) and r(t+x), which are generally 
known, and the predicted output of the ESO z x (t+x) and 
z 2 (t+x), which need to be obtained approximately. 

Prediction of f in ESO and RESO 

The primary method of predicting f is the Uor 2” rf -order 
55 Taylor Series approximation, i.e. 

Xt+lf)4l f )+ f( 64 or Equation (27) 
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2 Equation (28) 

5 

which requires the 1 st and/ or 2 nd order derivatives of f. To thi s 
end, a new form and a new way of implementation of 
extended state observer (ESO) 204 are introduced first, fol- 
lowed by the discussion on how to use them to obtain the 
predicted estimate off in Equation (27) and Equation (28). to 
Reduced-Order ESO 

In the observer-based control methods, the phase lag intro- 
duced by the observer decreases the phase margin of the 
control loop and is therefore undesired. To reduce the phase 
lag in the ESO 204, an embodiment of a reduced-order ESO 15 
(RESO) described here. 

Consider the plant 202 in Equation (9), if the derivative of 
the measured output y m (y m ), is given, a RESO can be con- 
structed with Zj and z 2 estimating x 2 =y and x 3 =f respectively. 
The correspondence between the augmented plant state space 20 
model and the RESO is shown as follows. 


18 

noise. In an alternative embodiment, a first order RESO can 
be constructed with the only state z tracking x 3 =f. The corre- 
spondence between the augmented plant state space model 
and the new RESO is shown as follows. 


0 1 0 
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0 0 1 

X + 
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u + 

0 

0 0 si 
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y = [1 0 0] x 

z = [ 0 ] z + [ :i; ]•<•! (y m - bu-z) 

A 

/ = z 

The shadowed blocks show how the gain matrices of the 
RESO are obtained from the augmented plant state space 
model. Equivalently, the RESO can be represented as 
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y = [1 0 0] * / = [0 1] i 


Equation (29) 
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The shadowed blocks show how the gain matrices of the 
RESO are obtained from the augmented plant state space 
model. Equivalently, the RESO can be represented as 

lt = A z z + B,u + L(y m - C z z) Equation (30) 

1 f = Vz 

where 


it = A zZ + B z u + L(y m -bu-C z z) Equation (35) 

1 f = Vz 

where A z =0,B z =0,L=l=co c) ,C z =l,V=l . The corresponding lin- 
40 ear PD control law 206 with this RESO is 

«=(g / ,^e,e,co c )+r- z )/6 Equation (36) 



where e=r-y m , e=r-y m . Note here again that the parameter- 
ized linear PD control law 206 can be replaced by the param- 
eterized NPD control law 506 to be introduced. 


C z =[l 0],V=[0 1], and 1 1 =2 co ( ,,1 2 =u) ( , 2 . 

With RESO replacing the ESO 204, the corresponding 
linear PD control law 206 in ADRC 400 is 

u= (gpd( e , e, <o t? )+r-z 2 )/ b Equation (31) 


A New ESO Implementation 

Consider the same second-order plant 202 as described in 
Equation (9). The ESO 204 described in U.S. patent applica- 
tion Ser. No. 10/351,664 can be divided into two parts. 


50 


where 


55 

e=r-y m ,e=r-z l Equation (32) 

If y m is clean it can be also fed back into the control law 
instead of z l , the Equation (33) changes to 


0 1 
0 0 


O'™ -zi) 


■/' 


Zi = \h{y m -Zi)dt 


Equation (37) 


e = r-y m , e = r-y m Equation (3 4) 

Note here that the parameterized linear PD control law 206 
can be replaced by the parameterized NPD control law 506, 
which will be introduced later. 

The order of ESO 204 is reduced further if y is also given, 
as the double derivative of y rn . According to U.S. patent 
application Ser. No. 10/351,664, the total disturbance is 
obtained as x 3 =f=y-bu and it is susceptible to measurement 


The first part is a Luenberger observer for the double inte- 
gral plant in Equation (13), and the second part shows that the 
estimation of the general disturbance can be obtained by 
integrating the observer error with a gain of 1 3 . This imple- 
mentation is equivalent to the ESO 204. 

When implemented in ZOH current estimator form, the 
Luenberger observer is 
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I z [k + 1] = 4>z[*] + Yu[k] + L p (y m [k] - z, [*]) 

\m = zlk] + L c (y m lk]-zdk]) 


Equation (38) 


ri t 

<d= , r = 

o 1 


(l+ft-2 Pi) /T 


2-2A 

(l+^ 2 -2A)/7’J’ 


l f l 3 s 3 4- {kpli 4- k d l 2 4- h)s 2 + Equation (44) 

1 ihph + k d l 3 )s 4- k\l 3 

C(5) = 

^ 5 3 + (/1 4- k d )s 2 4- (l 2 + k p + k d l\ — lfl 3 )s 

_ (s 3 + hs 2 + l 2 s + h)(s 2 + k d s 4- k p ) Equation (45) 

lfl 3 S 3 4- {kpl\ 4- k d l 2 4- l 3 )s 2 4- 
(&p/2 4- k d l 3 )s 4- &1/3 


where l x ,i 2 ,i 3 are observer gains defined in Equation (11) and 
k p and k d are the proportional -derivative (PD) gains defined in 
Equation (17). 

When the 2'“-order Taylor series approximation of f is 
1 5 used in control law 404, the equivalent 2dof transfer functions 


A = eT^cos (^w„r), h = ^ ° T ■ 


Similarly, the RESO in Equation (30) can be divided into 
two parts: 


zi =bu + h(y m -zi) 


Equation (39) 


: f l 2 (y m -i 


where the first part is a first order Luenberger observer, and 
the second part is an integrator. When implemented in ZOEI 
current estimator form, this Luenberger observer is 


Zi [k + 1] = 4>Zi [k] + r«[*] + L p (y m [k] - Zi [£]) Equation (40) 
Zi M = Zi [k] + L c (y m [k] - z, [*]) 


ijh 4 , i}hi,) 3 

~T S + w' 3 + — f 


Equation (46) 


1 (kpli 4- k d l 2 4- l 3 )s 4- (k p l 2 + k d l 3 )s 4- ki l 3 

CW = 1 3 I, 

s 3 + 4 -k d — Js 2 + 

( l 2 f hh\ 

^/ 2 +k p + k d li - lfl 3 — Js 

(5 s 4- l\S 2 4- l 2 S 4- l 3 )(s 2 4- k d s + k p ) 

Wh 1 4/1/3 ) 

{kpli +kdh + h)s 2 + [kpli + kdh)s + k\l 3 


Equation (47) 


Equivalent 2dof Transfer Function with Prediction in 
RESO 

The equivalent 2dof transfer functions of ADRC-RESO 
with prediction in f are as follows. 


where 0=1, T=bT, L=L p =l-$, and $=<r 2m ° T . 

Prediction of f 

With the new implementation of ESO 204, f=Jl 3 (y m -z 1 )dt, 

f=l 3 (y m -z 1 ) and f=l 3 (y m -z 2 ), and this enables us to determine 
f(t+ly) based on the D r or 2" -order Taylor series approxima- 45 
tion shown in 


Kt+lJSl&n-ztft+lJ&n-Zi) 


Equation (41) 


lfl 2 S? + (k p 4- k d l\ 4- l 2 )s 2 ■ 
1 (kpli + k d l 2 )s + k p l 2 
b s 2 + (k d 4- /i — l/l 2 )s 

(s 2 4- liS 4- l 2 )(kp 4- k d s 4- S 2 ) 
lfl 2 S 3 4- (kp 4- k d li + l 2 )s 2 4- 
(kpli 4- k d l 2 )s + k p l 2 


Equation (48) 


Equation (49) 


/(r + //)« k(y m -zi)dt+ifi 3 (y m -zi)+ -7r(y m -Z2) 


With RESO, f(t+L) is obtained as 


Xt+ijhfh (y m -zi)dt+i/ 2 (y m -zi) 


Equation (42) 


Equation (43) 


Equivalent 2dof Transfer Fimctions with Prediction 
in ESO 

Using prediction in ESO 204, with the U'-order Taylor 65 
series approximation of future value off used in control law 
404, the equivalent 2dof transfer functions are 


where l^lj are observer gains defined in Equation (35). 

State Predictor Method for Predictive Extended State 
Observer 

Another embodiment for fo rmin g a predictive state and 
disturbance observer 304 module is to use a state predictor 
method to extend the performance of a baseline ESO 204 to 
predict state estimates x seconds in the future. The basic state 
predictor method approach was published by T. Oguchi, FI. 
Nijmeijer, Prediction of Chaotic Behavior, IEEE Trans, on 
Circuits and Systems -I, Vol. 52, No. 1 1, pp. 2464-2472, 2005, 
which is hereby incorporated by reference. The state predic- 
tor method is combined in the embodiments presented herein 
with the regular ESO 204 of the form 


z =Az+Bu+L(y-y) 


Equation (50) 
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where 


y(t+-d=y m (t)+y m (tyz 


Equation (53) 
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Equation (51) 
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to form the predictive ESO is in the form of 

z(t)=Az(t)+Bu(t)+L(y(t)-y(t)) 

y(t)=Cz{t-x) Equation (52) 

Then y(t)=z 1 (t-x) will track y(t) while z^tj^ft+x) is its 
prediction x seconds aliead. The ideas is that since the 
observer design ensures that the observer error goes to zero, 
or becomes very small, if that error is defined as the difference 20 
between y(t) and Zj(t-x) and z^t-x) will approach y(t) and 
z x (t) will approach y(t+x). Therefore z x (t) can be used as a 
prediction of y(t) x second ahead. 

Building a Predictive State and Disturbance Observer Using 
an Output Prediction Mechanism 25 

Perhaps a simpler method to obtain z(t+x) is to leave the 
observer unchanged but replace its input y m (t) with the pre- 
dicted output y(t+x) as shown in FIG. 5. This configuration is 
used both to obtain the predicted error and its derivatives as 
well as the predicted disturbance estimation. Also depicted in 30 
FIG. 5 is a prediction module computer component or more 
generally a prediction module 508. The prediction module 
508 is adapted to accept the measured system output (y m (t)) 
and calculated estimated values of the system output (y(t+x)) 
from the present time to a time (x) in the future. 35 

A state and disturbance observer 504, shown in FIG. 5, 
accepts the predicted and current values of the plant output 
(y m ) and the control commands (u(t)) to estimate the future 
estimated states (z M (t+x)) where (n) is the number of esti- 
mated states plus an estimate of the generalized disturbance 40 
(f„). The state and disturbance observer 504 in one embodi- 
ment is an ESO 204. In a second embodiment, the state and 
disturbance observer 504 is a disturbance observer (DOB) 
structure. In yet another embodiment, the state and distur- 
bance observer 504 is an unknown input observer (UID). It is 45 
apparent that the combination of the state and disturbance 
observer 504 in combination with the prediction module 508 
provides a similar output and performance a similar function 
within the framework of an output predictive ADRC control- 
ler 500 as the predictive state and disturbance observer 304 of 50 
the predictive ADRC controller 300. In this manner, the vari- 
ous embodiments of the predictive ADRC controller 300 and 
output predictive ADRC controller 500 are correlated. 

Also shown in FIG. 5, a control law 506 is provided that 
uses the estimated future states of the plant 110 coupled with 55 
the desired trajectory (r^'^t+x)) to generate a current control 
input (u(t)). The embodiment of the control law 506 depicted 
in FIG. 5 incorporates the additive inverse of the estimated 
future values of the estimated disturbance (f„). 

Combine the Output Prediction with Regular ESO 60 

To compensate for the delay, y(t) is replaced with the 
approximately predicted output y(t+x) as the input to the ESO 
204. Again, the approximated prediction, y(t+x), is obtained 
from the Taylor Series approximation. The 1 '"-order Taylor 
series approximation and the 2" ,J -order Taylor series approxi- 65 
mation of y(t+x) are as shown in the following equations 
respectively. 


y(I+T) = y m (t) + y m (t)T + 


.V m W 


Equation (54) 


Then the predicted system output and the current control 
signal are used as the inputs to the ESO 204 as follows. 


| z(t + r) = A z z(r + r) + B z u(t) + Uy(r 4 r) - C z z(r + r)) Equation (55) 
l/0 + r) = Vz(r + r) 


which will provide both the predicted states and the distur- 
bance estimation in z(t+x). 

Equivalent 2dof Transfer Functions 

When the 1 "-order Taylor series approximation is used to 
obtain y(t+x), the equivalent 2dof transfer functions are 


(h 4 k d h + likp ).s _ + 

(1 + TS ) 

1 \ {kpk + hkd) s + kph ( 

b 5(5- + (/] + k d )5 + (h + k p + k d h)) 


H(s) = 


( 5 s +/15 2 +l2S + li)(kp +kdS+S 2 ) 
((3 + kdh +hk p )s 2 + 

(1 4 rs) 

(k p h + hkd)s + k p l} 


Equation (56) 


Equation (57) 


When the 2" rf -order Taylor Series approximation is used to 
obtain y(t+x), the equivalent 2dof transfer functions are 

/(/ 3 +M2+4*,>)5 2 + V t 1 A Equation (58) 

^ 1 1 (kph + hk d )s + kpl$ jl 2 j 

Equation (59) 

f (h + k-dh + hk p )s 2 + V r 2 ) 

[ikpi 2+ kk d )s + kpi 3 J( 1+ts+ t 52 J 

Combine the Output Prediction with RESO 

With the RESO configuration, the prediction is obtained as 

y(t+x)=y m ([)+y m (t)t Equation (60) 

and the new RESO is constructed as 


b 5(5- + (/] + k d )s + f/ 2 + k. p + k d li)} 
(s 3 + /j5 2 + hs 4 h)(k p 4 k d s 4 5 2 ) 


Z(t 4 T) = A z z(t 4 T) 4 B z u(t) 4 4 t) - C z z(t 4 r)) Equation (61) 

f(t+T) = Vz(t 4 r) 


where 



C_=[l 0],V=[0 1], and 1 1 =2oo 0 ,1 2 =<jo 0 2 . which again provides 
the predicted state and dismrbance estimation in z(t4x). 
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-continued 
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Equivalent 2 dof Transfer Functions 


With the above the D'-order Taylor Series approximation 
y(t+x) is used as the input to the RESO, the equivalent 2dof 
transfer functions are: 5 


g lv d(e, e, ai c ) = 


R sign(x), \x\ > R 
x, |x| < R 


where 



(k p + kjli + / 2 )s + 

(1 +Ti) 

C kpii + !ikp)s + k p l 2 
s(s + (l i + k d )) 


Equation (62) 


10 


H(s) = 


(5^ + lyS + l^)(kp 4" k^S 4" 5^) 

(k„ + k^ly + li)s 2 + 

(1 + TS) 

(k p l\ + l2kd)S + k p l2 


Equation (63) 


15 


Combine the Output Prediction with RESO 

With the RESO configuration, the prediction is obtained as 

y(t+x)=yjt)+y m (t)i Equation (64) 20 

and the new RESO is constructed as 


t(t + T) = A zZ (t+T) + B z u(t) + L(yU + T)-C z z(t + Tj) Equation (65) ^ 
f(t + T) = Vz(t + T) 

where 

30 


35 

C z =[l 0],V=[0 1], and l 1 =2co () ,l 2 =a> 0 2 . which again provides 
the predicted state and disturbance estimation in z(t+x). 

Equivalent 2 dof Transfer Functions 

40 

With the above the D r - 0 rder Taylor Series approximation 
y(t+x) is used as the input to the RESO, the equivalent 2dof 
transfer functions are 



signer) = 


l,x >0 
-1, x<0. 


Equation (68) shows that, within the region lxl<R and 
ly l<R, this parameterized NPD control law 506 is a> c 2 e+2co c e, 
which is the same as the parameterized linear PD control law 
206 in Equation (15); outside this region, however, the NPD 
control law 506 is nonlinear. This new NPD control law 506 
can be used to replace the linear PD control law 206 in the 
original ADRC 200 algorithm to achieve better performance. 
In one embodiment the NPD control law 506 is used in 
conjunction with a predicted state estimate z„(t+x). In a sec- 
ond embodiment, the NPD control law 506 is uses only the 
current state estimate (z„(t). 

A General Form of Reduced-Order Multiple Extended-state 
ESO for an n ,h -Order Plant 

Although the Predictive ADRC 300 has been presented 
with an ESO 204 or RESO with only one extended state, the 
same technique can be readily extended to a general n' ; '-order 
plant with an ESO 240 or RESO that has multiple extended 
states. To facilitate such development, we present here a gen- 
eral form of the reduced-order multiple extended states ESO 
204 for an n' ; '-order plant of the form 

y (n>= f(y,y , .... y <n ~ 1 ),d,t)+bu Equation (69) 

Let the n states be x 1 =y, x 2 =x 1 =y, . . . , , and 

the h extended states be x„ +1 =f, x„ +2 =x„ +1 =f, . . . , x„ +h = 
X rl+ fj_ | — U /— 1 ^ . The plant 110 in Equation (69) can now be 
expressed in the augmented form: 



(k p + k d l[ + 6 ).r + 

(1 + rs) 

(k p h + l 2 k d )s + k p l 2 
s(s + ( 6 +k d )) 


H(s) = 


(,F + /js + /,)(lp+ kj.s y s ~ ) 
(k p + k d li + h ).s" + 

(1 + TS) 

(kpli + l 2 k, i ).S +/,,,/; 


45 

Equation (66) 


Equation (67) 

x 


A Parameterized Nonlinear PD Control Law, Replacing the 55 
Linear PD Controller 

The parameterized linear PD control law 206 g pJ (e,e,oj c ) 
used above can be replaced by the following parameterized 
Nonlinear PD (NPD) control law 506 g l ^ d (e,e,to ( .): 

60 


V = co z e + io c e. 


( . VR(K + 8|vl) -R . , , 

I u c e + ^ 

[ oj c (oj c e + 2 e), 


\y\>R 

\y\<R 


Equation (68) 
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or equivalently 


( x = A x x + B x u + Ef 


0 


0 
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0 
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u + 
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0 
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or in state space form 


J = As + B,u + Hx„ + ,^ m+ 1 - c z z) 
f=Vz 


where 


to 


Equation (70) 15 
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Equation (71) 


where 
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if m=h,B z =[0] mxl ; if 
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35 


O' 


O' 
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b 

, v = 
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.0. 

mx 1 

.0. 


I- 0 kn+h)xi Thus, when m>h the nonzero elements of B z and V are in the 

(m-hy* and (m-h+l)**. column respectively. 

For tuning simplicity, the observer gains l l9 . . . , \ m are 
C =[1 0 . . . 0] lx(w+A) , and the nonzero element in B x is in the chosen as 
n'* row. 40 


Assume that y, y, . . . , are measured or otherwise 

available as y m , y m , .... y m ( ' n+> '~ m \ then the approximate 
states Xj, x 2 , . . . , x n+h _ m+1 are available (If m=h, the last state 
can be obtained from its definition: x M+1 =f=y„/" , -bu). Also 4 , 
assume i n> does not change rapidly and can be handled by the 
observer. The ESO 204 order can be reduced to m, where 
h<m<n+h. With x n+h _ m+1 and u as inputs, and with z 1 ,...,z m 
estimating x n+h _ m+1 , . . . , x n+h respectively, the general ESO 
204 form is 
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denotes the Binomial Coefficient 
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Also, to make the transfer function derivation easier, define 

60 



65 

With the state z m _^ +1 in the ESO 204 to estimate f, a control 
law, u=(u 0 -z m _ A+1 )/b, is applied to reduce the plant to a pure 
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n rf '-order cascade integral plant y ( ” )= u 0 , which can be easily 
controlled using the parameterized linear PD control law 206 
of the form 




0<i<n+l, and the middle term, 


an be easily and reject disturbance more effectively, either by control law 

trol law 206 modification or ESO 204 modification. With a part of the 

plant dynamics known, the known part of the generalized 
disturbance can be generated and rejected in the control law 
5 directly, while the rest of the generalized dishirbance is esti- 
Equation (72) mated by ESO 204 and then rejected in the linear PD control 
law 206 as previously shown in the original ADRC 200 
framework. The other way to utilize the partially known plant 
dynamics is adapting the dynamics into ESO 204, where ESO 
to 204 still estimates the generalized disturbance. 

Consider the second-order plant 202 in Equation (8), with 
the knowledge of the plant information, part of the general- 
ized disturbance is known, denoted as f„(y,y). The rest is 
unknown part of the generalized disturbance, denoted as f„(y, 
15 y,d,t). The plant in Equation (9) can now be described as 

y=f n (y,y, d,t)+bu Equation (76) 

and there are two ways to take advantage of new information 
in f„(y,y) and incorporate it into ADRC 200. One way is to 
20 improve ESO 204 by combining f„(x) with A v x+B x u in the 
plant space description in Equation (76). The augmented state 
space form of the plant is 


is corresponding to the order reduction (it disappears when 
the inputs of ESO 204 are y m and u). In the case where the 25 
measured output y m is clean and its derivatives can be 
obtained from y m , ADRC 200 can also be applied as 


uo=Yj k ' r ° ” _ E ^ l) - Zm ~ h * 1 

;=i i=i 


72+ 1 

j =_ 1 

n+h-m 

n, m z kjs(j-» + 
j = 1 

m-h+l 

E l , .rn+h-l-m+jXl)m 

•'•n+h-m+ j* T J,m 

j i=l 


Equation (74) 


( x = A x x + Ef n {x ) 4 - B x u -\ 

\ y = c x x 


Equation (77) 


and the corresponding ESO 204 is 


Equation (73) 30 


The corresponding transfer function representation of 
ADRC 400, in a 2dof structure (as shown in FIG. 4) is 35 


t = A x z + Ef n (z ) + B x u + Uy m - C x z) 
f=Vz 


Equation (78) 


Then a new ESO 204 can be constructed to estimate the 
generalized dishirbance f based on the new state space 
description of the plant 110 in Equation (77). For example, 
consider the plant y=-a 1 y+bu where a : is given, i.e. 
f„(y,y)=- a iy. The known part of the generalized disturbance 
is obtained from the model information: f„=-a 1 y=-a 1 (-a 1 y+ 
bu). In this case where f„(y,y) is a linear function, the aug- 
mented state space representation of the plant 110 can be 
written as 


n+h-m 

r„, y k jS v-» - 

j = i 


Equation (75) 


x = A x x + B x u + Ef u 
y = c x x. 


Equation (79) 


1 J= l 

G c (s) = 7 — — 

b m t +l 


J,m 
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where A x = 

0 0 1 

, B x = 
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U Vb"'- 


55 and the new ESO is z=A x z+Bm+L(y m -C x z), where B z =0 

b-a 1 b] r 

The other way of taking advantage of the additional infor- 
mation in f„(y,y) is to use it to cancel known plant dynamics 
in the control law directly. There are three ways to do this: 


0<x<y. 

ADRC with Partially Known Plant Dynamics 

ADRC controller does not require the detail plant model 
information. Flowever, if the whole or part of the plant 65 
dynamics is known, various embodiments of the ADRC con- 
troller can be modified to utilize the information to estimate 


M=(g(e,e,(O c )+r-/ J (y 77 J , i )> 77 J )-Z 3 )/^ Equation (80) 

or y rn and y m can be replaced by z x and z 2 in the control law, 
u ={g{e,e^ c ) + r-f n {z h z 2 )-z 3 )lb Equation (81) 

or y m and y m can be replaced by r and r respectively and 
rejected in the control law in a feedforward manner. 


«=fe( e,e,(i> „)+f-f„(r,r)-z } )lb 


Equation (82) 
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where g(e,e,o) c ) is the parameterized PD controller 206 or 4. The controller of claim 2, wherein the control module 

NPD controller defined in Equation (15) or Equation (68) comprises a non-linear control law given by an equation: 
respectively. Of course, any combination of the above three 
methods are available for use by one of ordinary skill in the 

art. 5 y - oj^e 4- to c e. 


INDUSTRIAL APPLICATIONS 

The present active disturbance rejection controller, and its 
associated processes, methods and devices as disclosed 
herein possesses a number of unique attributes and industrial 
applicability, including for example, utility as a controller for 
tracking control, web processing applications, and jet engine 
control applications as described in greater detail above. 

For example, in the case of tracking or motion control 
applications, the various predictive methods presented 
herein, particularly the ones that use the knowledge of the 
plant dynamics, offer reduced phase lag in the observer and 
compensate for the time delay in the plant itself, all of which 
especially benefit tracking applications that employ a desired 
output trajectory, sometimes called motion profile. 

CONCLUSION 

While various embodiments of the present system and 
method for feedback control of systems have been described 
above, it should be understood that the embodiments have 
been presented by the way of example only, and not limita- 
tion. It will be understood by those skilled in the art that 
various changes in form and details may be made therein 
without departing from the spirit and scope of the invention as 
defined. Thus, the breadth and scope of the present invention 
should not be limited by any of the above described exem- 
plary embodiments. 

What is claimed is: 

1. A controller, comprising: 

a predictive state and disturbance observer module config- 
ured to receive a sensor signal and a disturbance adjusted 
control signal, and to output extended state estimate data 
representing a state of a plant, an extended state of a 
plant system dynamic and an external disturbance of the 
plant, and a predicted extended state estimate, wherein 
the sensor signal is based on an output associated with 
the plant; and 

a control module configured to output a control signal 
based on the extended state estimate data, the predicted 
extended state estimate, a trajectory and a trajectory 
prediction; and 

wherein the extended state of the plant system dynamic is 
adapted to cause an input-output characteristic of the 
output and an associated input of the plant to behave as 
a double-integral plant with a scaling factor. 

2. The controller of claim 1, wherein the predictive state 
and disturbance observer module comprises a system output 
predictor and an extended state observer. 

3. The controller of claim 2, wherein the system output 
predictor is configured to predict a future value of the sensor 
signal using a Taylor series approximation. 


( . V R(R + 8|y|) - R 

x= \ co c e+ sign(_y), \y\ > R 

\ oj c {(jj c e + 2e), |;y| < R 

( R sign(*), \x\ > R 
g^e,e,a> c ) = I 


( 1, -v > 0 
signU) = 1 


and where u) c is a frequency of the controller, e is an error, 
R is a maximum control signal, and % ndp is the control 
signal. 

5. The controller of claim 2, wherein the sensor signal 
represents at least a first derivative of the output of the plant, 
and wherein the extended state observer is a reduced order 
extended state observer. 

6. The controller of claim 1, wherein the predictive state 
and disturbance observer module comprises a predictive 
extended state observer. 

7. The controller of claim 1, wherein the predictive state 
and disturbance observer module comprises an observer 
model of a dynamic of the plant. 

8. The controller of claim 1, wherein the control module 
further comprises an additive inverse model of a function (f„) 
that comprises an estimate of a dynamic of the plant and the 
external disturbance. 

9. A computer-implemented method for controlling a 
plant, comprising: 

estimating an extended state for the plant using a sensed 
output of the plant, wherein the extended state comprises 
a state of the plant and a total disturbance of the plant, 
and wherein the extended state is adapted to cause an 
input-output characteristic of the plant to behave as a 
double-integral plant with a scaling factor; 
predicting a change in the extended state to generate a state 
prediction and an extended state prediction; 
specifying a trajectory for an output of the plant to follow, 
wherein the trajectory comprises a desired trajectory and 
a prediction of a future desired trajectory; 
applying a control law to the desired trajectory and the 
prediction of the future desired trajectory, the extended 
state, and the extended state; and 
generating a control output based on application of the 
control law. 
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